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Abstract. We study the possibility that charge breaking minima occur in the Zee model. We 
reach very different conclusions from those attained in simpler, two Higgs doublet models, and 
the reason for this is traced back to the existence of cubic terms in the potential. A scan of 
the Zee model's parameter space shows that CB is restricted to a narrow region of values of the 
parameters. 



1 Introduction 

In models with several scalar doublets there is the possibility that minima which break charge 
and/or colour symmetries occur. That happens when scalar fields which are not singlets of 
charge and/or colour acquire a non-zero vacuum expectation value. If those minima are deeper 
than the standard ones, which respect the SU{3)c x U{l)em gauge symmetries, the theory 
might tunnel to them, and find itself in a vacuum that breaks charge and/or colour. As such, 
the combinations of parameters of the theory that give rise to those minima can be excluded, 
and the predictive power of the theory is increased. These charge and colour breaking bounds 
were first introduced by Frere et al and have been extensively used in supersymmetric models, 
to impose limits on the masses of the supersymmetric partners, for instance [2]. Recently |31 the 
same idea was applied to two Higgs doublet models (2IIDM) without explicit CP violation and a 
very interesting result was deduced: if a minimum that respects electric charge conservation and 
also the CP symmetry (we will call this minimum the "normal" one, from this point onwards) 
exists, then it is the global minimum of the theory. In fact, it was shown that only charge 
breaking saddle points can occur in that case, and the value of the potential in those points is 
always superior to its value at the normal minimum. In this paper we will consider the simplest 
generalisation of the 2IIDM, in which we add an SU{2) singlet to the theory. This is the Zee 
model which was first conceived as a means to explain the very low masses of the neutrinos. 
The Zee neutrino mass matrix structure can be shown to arise in a more general context, through 
mixing of dimension five operators ^ . We will show that in the Zee model charge breaking (CB) 
minima can and do develop in parallel with normal ones, and sometimes those CB minima are 
the deepest ones. Then, unlike what happens for the 2HDM, tunneling to CB minima is actually 
possible even if the theory starts up in a normal minimum. We will investigate the values of 
the potential's parameters for which CB occurs, and show that CB is mostly circumscribed to 
a specific region of the parameter space. 
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2 The Zee model 



The Zee model jl] has the great advantage of providing a natural reason for the smallness of the 
neutrino masses. Majorana mass terms for the neutrinos are generated as radiative corrections 
at the one-loop level. In addition, the model can easily be extended ^ to incorporate a sterile 
neutrino, in such a way that this four-neutrino version could also accommodate the LSND data. 
Unlike the seesaw mechanism, the Zee model does not require new physics at an extremely 
large scale. But it does require a Higgs sector larger than the SM's, i.e., it requires two scalar 
doublets, $a (a=l, 2) with hypercharge Y = 1 and a charged singlet field $o with hypercharge 
Y = —2 ^. In this work we are interested in determining the possible vacua of the theory, i.e. we 
want to study the scalar potential of the model. The 2HDM scalar potential is given by a sum 
of terms quadratic and quartic in the fields and ^2- There is an extensive literature about 
this potential (see, for instance, so we simply summarise some of its features: (a) the most 
general SU{2) x U{1) gauge invariant potential, which explicitly breaks CP conservation, has 14 
independent real parameters, that can be reduced to 11 with a judicious choice of basis jH]; (b) 
requiring that the 2IIDM does not explicitly break CP one is left with a 10 parameter potential 
(though with a choice of basis we can reduce this number to 9 independent parameters); this 
is the potential studied in reference Pj, end it has flavour changing neutral currents (FCNC); 
(c) if one wishes to prevent the occurrence of FCNC's - or indeed the possibility of spontaneous 
breaking of charge or CP symmetries -, it was shown in reference [2] that the potential has to be 
invariant under a Z2 or U{1) symmetries; one then obtains two different 7-parameter potentials. 

The Zee model usually considered in the literature is based on the 2IIDM described in point 
(c) above. A priori this seems like a good idea, since in this manner we prevent potentially 
dangerous FCNC's and limit the number of free parameters of the theory. The problem is that 
recent analysis |in[ lllj of experimental results in solar and atmospheric neutrinos, namely from 
the SNO experiment, seems to indicate that the minimal Zee model cannot account for the data. 
It has been suggested that extensions of the Zee model with extra scalars might do a better job 
of describing the neutrino results (for instance, see jlH I12j). It has also been shown that a 
Zee model based on the 2HDM potential described in point (b) above could accurately reproduce 
the experiments, and that is the case we will study here. An obvious problem with this approach 
is the fact that we have increased the number of free parameters in the theory. So, any tool that 
might be useful in reducing the parameter space of the model will be welcome. For instance, 
the authors of refs. ^1] derived mass bounds on the Higgs boson using the triviality condition 
and (standard) vacuum stability, both for the Zee and two higgs dublet models. If we can show 
that for certain values of the parameters deeper CB minima occur, we will be able to impose 
limits on those parameters and improve our understanding of the theory. To achieve this, we 
must look at the minimisation conditions for the normal and CB minima. 

The fields {^q , , ^2} have ten real component fields, which we number as 



This convention might seem odd, but it is very convenient for writing the mass matrices of the 
scalar particles of the theory. We can build five SU{2) x f7(l) quadratic invariants with these 
fields, which are: 




(1) 



X2 = \^2? =^1+^1+^1+^1 



Our convention for the electric charge operator is Q = Ts + Y/2. 



X4 = lm{^\^2) = VlfA: - f2V3 + V5f8 - f6V7 ■ (2) 

The construction of the scalar potential is thus very simple, one must simply include all terms 
linear and quadratic in the variables Xj. We define the charge conjugated states as = 
$*. Under this transformation X4 — >■ — X4 and all the remaining Xj are invariant, hence a CP 
conserving potential is 

Vm = aoxo + ai xi + 02X2 + 03X3 + 

600 + 611 xl + 622 + 633 X3 + 644 xl + 

601 xqXi + 602 X0X2 + 603 X0X3 + 612 X1X2 + bi3 X1X3 + 623 a;2a;3 • (3) 

The ai parameters have dimensions of mass squared, the bij are dimensionless and all the pa- 
rameters are real. For future convenience let us introduce here a matrix notation similar to the 
one used in ref. [3]. Let X be the vector with components (xq , xi , X2 , X3 , X4) and A and B 
the vector and symmetric matrix given by 
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Therefore, the potential Vm is given by 

Vm = A^X + \x'^ BX . (5) 

However, there is also an SU{2) x U{1) invariant cubic term that gives an extra contribution to 
the potential, namely 

Vc = M ($f icj2 ^2) ^0 + h-c. 

= 2M [{^PUPq - ip2'~P8 - ^3V>5 + ^i^r) V>9 - {^1^8 + V2V6 - ^3^7 " ^i^5) ^lo] , (6) 

where M is a real coupling with dimensions of mass. The scalar potential of the Zee model is 
thus V = Vm + Vc- 

The "normal" minimum, which preserves charge and colour conservation, occurs when the 
fields 995 and ipQ have non-zero vevs, vi and V2 respectively, such that vf + V2 = (246 GeV)^. At 
this minimum, then, the vector X takes the values X^ = (0 , , f | , V1V2 , 0), and the value 
of the potential is given by jS] 

Vn = Xn = ^ {aivj + a2V^ + a3ViV2) < . (7) 

The inclusion of the charged field introduces some changes in the theory's scalar mass 
spectrum. If one computes the second derivatives of the Zee potential with respect to the fields 
ipi at the normal minimum, we see that the neutral sector of this model is identical to that of 
the 2HDM. The reason for this is quite simple: $0 is an electrically charged field and so there 
are no mixing terms with the neutral components of $1 and $2- A different thing happens in 
the charged sector, where there is a mixing among the charged fields caused by the cubic terms 
Vc- The charged higgs mass matrix is now given by (from the left to the right, the columns of 



this matrix correspond to the fields ipi, ip2, ^3, (fi, ^9 and <y9io) 
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where we have introduced the notation Vl 
normal minimum and one has 

V2 



^/ = - — Vi 

^ 2vi ^ 



Vi 



dVM/dxi. The Vl, of course, are evaluated at the 

(9) 
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The matrix (jSJ is clearly a generalization of the analogous matrix obtained in reference ^ for the 
2HDM. This matrix has two zero eigenvalues (corresponding to the charged Goldstone bosons) 
and two doubly degenerate non-zero ones, given by 



^1,2 



Vi + Vl + Vi ± J {VI + Vi - Vif + 4M^{vl + 



(10) 



3 Charge breaking minima 

For charge breaking to occur, some charged scalar field has to acquire a non-vanishing vev. One 
might think that the case {ips ^ , (pi = 0, Vj^a} (or others, equivalent to this one, where a 
single charged field of the doublet has a non-zero vev) would give rise to a CB minimum. In 
fact, it is simple to show that such a choice of vevs for any SU{2) doublets does not actually 
break charge conservation jl5j . However, giving a vev to the singlet definitively gives the 
photon a mass. Let (/jg = 7 7^ (or ipiQ 7^ 0, which is equivalent) and all the other fields zero. 
In this case only xq = 7^ is different from zero, and a trivial minimisation of the Zee potential 
shows that 7^ = — oq /2 6oo- For the potential to be bounded from below we must have 5oo > 
and so this minimum only exists if oq < 0. In that case the value of the potential is 

Then, charge breaking does not occur if Vcbi > Vjy, which implies that the theory's parameters 
must obey the inequality 

^00 > -I r- 4- > • (12) 

2 aivf + a2 V2 + as vi V2 

Given {ap , ai , 02 , 03 , f 1 , ^2}, we have thus a restriction on how small the parameter 600 may 
be. This CB minimum is rather trivial, and it seems to impose constraints on only a few of 
the theory's 16 free parameters. In ref. !3. we considered a CB minimum of the form Lp^ = v'l, 
ipQ = v'2 and 933 = Q. It is easy to see, though, that for the most general choice of parameters 
of the potential this CB minimum is not always possible. The minimisation conditions are 
dV/difi = 0, Vj=i^...^io and we wish to find if a solution with only {ip^ , ip^ , ipg} 7^ is possible. 
From dV/d(pg = we obtain 2Mv[a = 0, which, for a 7^ 0, implies either M = (a case of 



A 



limited interest, as the model almost reduces to the 2HDM) or v'l = 0. From dV/dip^^ = 0, 
with if^ = v[ = and tpQ = V2 0, we obtain 



{ , (13) 

[ 2a2v'2 + 4622^2 = 

which can only be satisfied for particular combinations of parameters, namely 03/623 = 02/2 622) 
with a2 < 0. 

Let us study the general case and take the parameters of the theory as being uncorrelated. 
The conclusions we will reach should therefore encompass all the special cases, such as the 
examples we have just discussed. A simple investigation of the minimisation conditions quickly 
tells us that the simplest CB choice of vevs, besides the trivial one considered above, is ipi = p, 
ip3 = a, ip5 = v[, (fQ = v'2 and (/?9 = 7. This gives xq = xi = v'^ + p^, £2 = v'^ + a^, 
X3 = Viv'2 + pa and £4 = (we use the hat to distinguish these Xi from their counterparts at 
the normal minimum). In reference |S] a very useful relation between the values of the 2HDM 
potential at the normal and CB stationary points was derived, namely 

VcB - Vn = (-^) [M^2 - V2vif + a'vl] 

= 2{vlXvl) + a-vl] , (14) 

where {v'l , ^2 , a} are the CB vevs and = dV/dx^ is evaluated at the normal stationary 
point. Mjj± is the squared mass of the charged scalars at the normal stationary point. The 
second line in the equation above makes it plain that if the normal minimum exists then we will 
necessarily have Vcb — Vn > 0, and tunneling to a deeper CB minimum is impossible. Let 
us now generalise this result for the present case. As a first step we define the vector V' , with 
components V- = dVm/dxi evaluated at the normal minimum, i.e., 

V' = A + BXn . (15) 

At the CB stationary point we define the vectors Xqb, with components XcBi = Xi, and 
V' = V- = dVM/dxi\^^, that is, 

V' = A + BXcB ■ (16) 
It is easy to see, from the CB stationarity conditions, that we have the following relations: 

Vi = . VI = . Vi = -^i'i . VI = lilf^ia^J , V-i = . (17) 

2xq 2x3 2x3 Vc 

Vc is now evaluated at the CB stationary point and is equal to Vc" = 2M7 (pwg — av'^). Since 
the potential is a polynomial with quadratic, cubic and quartic terms, it is trivial to show that 
its value at any stationary point is given by half of the quadratic terms plus a quarter of the 
cubic ones. Then, it follows that 

Vcb = \xIj,A + -^Vc . (18) 
From eqs. ()15p and ()16|) we obtain 

r xi^v = Xl^A + XI^BXn 
\ XlV' = XlA + XlBXcB 



Because the matrix B is symmetric, we have X^^ B Xj^ = Xjj B Xqb- Then, subtracting 
these equations, and using equations Q and ((TH|) we obtain 



xIbV - XlV' = xIbA - XlA = 2{VcB - Vn) - \vc . (20) 
On the other hand, from eqs. Q, we have 

X^sV'= + + P^)V; + + a-')V^ + {v[v'^ + pa)y3' 

= 7'^o - ^ [iv'iV2 - v'^vi)^ + {vip - V2a)^] , (21) 

and (fTTj) give 

Xj^ V' = vl Vi + vl Vi + viV2 Vi 

= - [M^2 - v'^vif + {vip - V2af] . (22) 
2x3 

Replacing these expressions in (|2U|) and using (|17j) . we finahy obtain 

2(l/cB - Vn) = 7' - V^o) + M"^ - — ) [(^'1^2 - v'^vif + {vip - V2af] . 

V / z \ X3 a;3 / 

(23) 

There is a striking similarity between this expression and the one found for the 2HDM potential, 
equation (fTl]) . However, whereas ((TH) establishes unequivocally that Vcb — Vn > 0, for the Zee 
model this is not possible. Depending on the values of the parameters of the theory, we may well 
have Vq > Vq and V^/x^ > V^/x^, or exactly the opposite. Hence, for certain combinations of 
parameters, it is a priori possible that a CB minimum deeper than the normal minimum exists. 
It is the existence of the cubic terms in the potential that makes the simultaneous existence of 
normal and CB minima possible. In fact, if M = eqs. (|17|) show that = Vq = 0, and 
that V^' > and —V^/x^ > are the values of the two charged higgs masses at the normal 
minimum. So we trivially obtain, from eq. (|23|) . Vcb > Vn. Notice that the presence of the 
cubic term makes the Zee model more similar to supersymmetric theories than to the 2HDM, 
and one knows that CB can occur in SUSY theories. 



4 Numerical results 



In view of the result that we have obtained in the previous section, let us try to determine what 
combinations of parameters {M , Oj , bij} could give rise to CB minima deeper than the normal 
one. It is impossible to solve analytically the CB minimisation conditions. Even numerically, 
this turns out to be a difficult task. So, we adopt the following strategy: we choose the values 
of the parameters such that V has a normal minimum. For the same parameters we now take 
F as a function of {v'l , V2 , p , a , j}. Starting from an initial guess for the values of these CB 
vevs, we numerically sample the value of the potential in its neighborhood and continue along 
the direction for which V decreases. At the end of this evaluation, we compare the value of V 
that we have reached with Vn. If it is smaller than Vn, then this set of parameters is rejected, 
because it could give rise to CB. We tried to overcome the dependence of this method on the 
initial guess by performing multiple runs of the algorithm with different initial conditions. 

The potential depends on 16 independent parameters. We randomly generate values for 
14 of these parameters, leaving oi and a2 out of this procedure. The value of tan/3 = V2/V1 



is also used as input (with 0.5 < tan/3 < 100). We impose some restrictions, though. The b 
parameters are chosen as being of the same order of magnitude, since it would not seem natural to 
do otherwise. Here we were inspired by SUSY models, in which all bij are of similar magnitudes. 
Likewise, we take the {M , a,} parameters to be of the same order. The intervals of variation 
for the parameters are: -3 < M < 3 TeV, -1 < Oj < 1 (TeV)^ and 0.01 < bij < 10. Having 
chosen these parameters, we use the values of f i = v cos (3 and V2 = v sin (3 (with v = 246 
GeV) and the minimisation conditions of the normal minimum, eqs. Q, to determine the values 
of fli and 02, given by 

I 

ai= - - — Us V2 + 4: bu vf + 2 (612 + 633) "Ui + 3 613 vj V2 + 623 ■^2] 

a2= - [03 + 4 622 + 2 (612 + 633) vlv2 + 3 623 vi + 613 vf] . (24) 

2V2 

With the values of these parameters we compute the masses of the scalar particles, using the 
expression ^ for the charged scalars and the formulae (37) and (38) of reference |3] for the 
neutral scalars. We accept only those sets of parameters for which (80)^ < m^, < (1000)^ 
(GeV)^. This interval of variation for the squared masses was selected so that they are positive, 
thus ensuring we are at the normal minimum, and larger than the values already ruled out by 
experimental searches and inferior to the usual triviality bounds. 

It is possible that, for some combinations of bij, the quartic terms of the Zee potential 
become arbitrarily negative when some of the fields tend to infinity. This corresponds to a 
potential which is unbounded from below (UFB) and this combination of bij should be rejected 
on physical grounds. In ref. |3] we gave some necessary conditions on the 6's to prevent the 
occurrence of UFB directions in the 2HDM, and similar conditions were deduced and used in 
this work. The procedure to find them is straightforward: take one or more of the (fi fields 
to infinity and determine for which values of bij we might obtain V —00. The use of 
polar /spherical coordinates is very helpful to study these limits. The most trivial conditions 
one obtains in this fashion are {600 j ^11 i ^22} ^ 0. The most stringent UFB bounds were 
found taking {ipi , ip^ , 999} — s- 00 and parameterizing these fields as tpi = r sin a cos 6*, = 
r sin a sinO, ipg = r cos a, with r +00 and {a, 6} arbitrary angles. Requiring that the 
limit of the quartic terms of the potential is larger or equal to zero produces the condition 
fia,9) < g{a,6), M^fi, with 



f[a,e) = 



^ sin^ 2a + (613 cos^ 9 + 623 sin^ 6) 



2 

g{a, 9) = 600 cos^ a + 
1 



611 cos^ 9 + 622 sin^ ^ + ^ (^12 + 633] sin^ 29 ] sin^ a + 



^ (601 cos^ 9 + bo2 sin^ 9) sin^ 2a . (25) 

We checked, for several values of {a , 9}, if the values of the b's obeyed the previous inequality. 
The UFB conditions found in ref. |2i for the 2HDM are still valid for the Zee model and have 
been used. Furthermore, it is easy to obtain two additional conditions. 



boi > -2 V^oob 
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bo2 > -2 Vboob22 . (26) 

Although these are only necessary conditions, we found that they eliminate most of the UFB 
potentials. Nevertheless, we will see that UFB potentials still appear in our analysis, but only 
for a small subset of parameters. 



Having required these conditions on the b parameters and rejected those sets of {M , Oj , bij} 
that did not produce scalar masses within the given interval, we built a set containing nearly 
103000 different combinations of parameters. Although this a small part of the total parameter 
space, we hope that it will serve to illustrate the general CB properties of the Zee model potential. 
For each "point" in this parameter space, we computed the value of the potential at the normal 
minimum, Vn, numerically obtained the value of Vcb and compared both these values. Despite 
our effort to enforce the UFB conditions, we have obtained UFB potentials for roughly 2% of 
the points. 

CB was found to occur, mostly limited to a specific region of parameters. In fact, the 
majority of CB points occurred for large absolute values of the M parameter {\M\ > 500 GeV), 
which is easy to understand, given our earlier conclusion that the existence of the cubic terms in 
the potential is intimately linked to the possibility of CB. The region where CB was most likely 
to occur has lower values of tan/3 (inferior to about 10), negative values of the parameters 633, 
644 and as (this last one always superior to —(800 GeV)^) and positive values for ai and bi2- 
Having varied the b parameters in an interval spanning three orders of magnitude (from 0.01 
to 10) we found that CB occurs less frequently for the larger values of the Vs. Although these 
observations are interesting in themselves, looking for the possible impact of CB on physical 
parameters, such as scalar masses, is of much greater interest. We observed that CB occurred 
for comparatively smaller values of the lightest neutral (CP even) higgs scalar, mostly below ~ 
350 Gev with the largest concentration occurring for 90 < < 130 GeV. This is illustrated 
in figure where we plot the values of Mh versus Mh- The area delimited by the dashed lines 
shows the entire zone filled out by our ~ 10^ sample of the parameter space. Those points for 
which CB occurs are marked with a cross. It is easy to see that CB occurs mostly for a narrow 
region of values of Mh. Similar "trends" in the values of the masses are found for the other 
scalars of the theory. CB is concentrated on the relatively small values of the lightest charged 
higgs (M^± < 300 GeV), meaning larger values of the heaviest charged scalar {M^± > 500 
GeV), and for values of Ma, the pseudoscalar mass, usually larger than 200 GeV. 

This initial foray into the CB analysis of the Zee model parameter space allowed us to 
identify the areas for which it was more likely to occur. Our next step was to generate a second 
set of parameters (about 53000, this time), following the same philosophy but with some extra 
restrictions: namely, tan (3 < 10, 0.01 < bij < 1 and imposing M^ < 500 GeV. Our intent 
was to narrow the search for CB to those areas where it is more likely to occur, to gauge its 
importance in terms of percentage of parameter space rejected. We found that CB occurred 
for about 5% of this second parameter space, and confirmed the "trends" found earlier for the 
parameters and masses of the theory. As an illustration of this, the reader can see in figure[21how 
CB occurs for smaller values of the lightest charged scalar (below 300 GeV, as mentioned above) 
and mostly for higher values of the pseudoscalar mass. Again, in this figure the area delimited 
by the dashed lines shows the entire range of allowed values. Unfortunately, these plots do not 
mean that, after having eliminated all the points of the parameter space for which CB occurs, 
we would be left with "holes" in this space. If it were so, we could make definite restrictions, for 
instance on the higgs' masses. But this is not the case: given the size of the parameter space, 
there are many combinations of {M , aj , bij} which produce the same scalar masses. In other 
words, if in figure ^ we showed that CB occurred for 90 < Mh < 130 GeV, that does not 
mean that all Zee model points which produce masses in that interval were eliminated. On the 
contrary, there are many left which produce such masses. This means that we cannot produce 
a definite statement about the areas of parameter space for which CB occurs. We can, however, 
say with a good degree of certainty where CB does not occur: the regions with high values for 
the masses of the lightest neutral and charged higgses {Mh > 500 GeV, M^± > 300 GeV) are 
almost entirely "clean" of CB. 



5 Conclusions 



We have shown that charge breaking may occur in the Zee model, and have developed here the 
algorithm for finding the regions of parameter space where that may happen. The values of 
the lightest neutral higgs mass for which CB is more likely to occur are specially interesting, 
given that they correspond to the values that are currently being probed by our accelerators. 
Our analysis was very general, considering large regions of parameter space with randomly 
generated parameters. As was mentioned in the text, the minimal Zee model does not seem 
to be compatible with current neutrino data. This work focused on more general forms for the 
two doublets plus one singlet scalar potential, but extensions of the Zee model, with a larger 
scalar sector, are also of interest. The possibility of charge breaking would also arise in those 
models, and it would be interesting to perform CB studies on those theories. The results we 
obtained here cannot, clearly, be generalised to those other models, and in fact there is the 
possibility that, given a particular scalar content, the conclusions reached will be very different. 
For instance, in ref. jl6j a model with a hyperchargeless triplet majoron was studied, and it was 
shown that the vacuum of that theory always broke charge conservation. Extra scalar content 
of a more "usual" nature - singlets and dublets, charged or not - would make the Zee potential 
even more similar to the super symmetric one, and thus we would expect that charge breaking 
would still occur, as it does in the SUSY case. 

Given the results encountered for our study of the non-minimal Zee model, we feel that 
further studies in this area are warranted: specific embodiments of the Zee model, with values 
of parameters specially tuned to explain neutrino observational data, should be tested to verify if 
they are safe from CB. This is not usually done, since in most applications of the Zee model |in| 
[TTl [T^ one does not consider in detail the scalar sector of the theory beyond its spectrum of 
masses and mixing angles. However, such a study might help narrow down the range of variation 
of the theory's parameters and, perhaps, shed some light on the generation of neutrino masses. 
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Figure 1: Plot of the mass of the lightest CP-even higgs neutral scalar versus the mass of the 
heaviest; the area delimited by the dashed lines shows the whole region of the parameter space. 
The crosses denote those "points" for which CB occurs. Most of these are circunscribed to a 
narrow interval of low values of M/j. 
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Figure 2: Plot of the mass of the lightest charged higgs scalar versus the mass of the pseudoscalar; 
the area delimited by the dashed lines represents the whole of the parameter space, and the 
crosses those points for which CB occurs. 
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